Anomalous fluctuations of the condensate in interacting Bose gases 
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We find that the fluctuations of the condensate in a weakly 
interacting Bose gas confined in a box of volume V follow the 
law {5 No) ~ V 4 ^ 3 . This anomalous behaviour arises from the 
occurrence of infrared divergencies due to phonon excitations 
and holds also for strongly correlated Bose superfluids. The 
analysis is extended to an interacting Bose gas confined in a 
harmonic trap where the fluctuations are found to exhibit a 
similar anomaly. 
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The recent experimental achievement |lj of Bosc- 
Einstein condensation (BEC) in dilute atomic vapours 
has stimulated a great revival of interest in the theoret- 
ical study of this phenomenon (for un update list of ref- 
erences see Q). Among the several intriguing questions 
associated with Bose-Einstein condensation the problem 
of fluctuations plays an important role in the more gen- 
eral context of statistical mechanics (see for example Ref. 

The fluctuations of the condensate have been studied 
in a systematic way in the case of the ideal Bose gas 
in a box and, more recently, in the presence of 

a harmonic trap In the ideal Bose gas one finds 

a dramatically different behaviour between the canon- 
ical and the grand canonical ensembles. In particular 
the results of Ref. ||,[| show that in the canonical en- 
semble the fluctuations of the number of particles in the 
ground state of a box of volume V behave as V 4 I Z . This 
law ensures that the relative fluctuations ((SNq) 2 } /V 2 of 
the condensate vanish when V — > oo, differently from 
what happens in the grand canonical ensemble, where 
((SN Q ) 2 ) = N$ + N„ ~ V 2 . The results obtained in 
the canonical ensemble, however, reveal that the fluctu- 
ations of the condensate are not linear in V, in contrast 
to the fluctuations of most physical observables. One 
might think that this anomalous behaviour is a pecu- 
liarity of the ideal Bose gas and that, as soon as one 
includes interaction effects, the fluctuations of No take a 
normal behaviour, linear in V. In this letter we address 
explicitly this question and we find that, despite the fi- 
nite compressibility of the system, the fluctuations of the 
condensate still behave as V 4 ^ 3 . 

Let us start our investigation in the framework of Bo- 



goliubov theory of a uniform weakly interacting Bose gas 
confined in a cubic box of size L. According to Bogoli- 
ubov theory Q the particle operator can be expressed 
in terms of the quasi-particle creation and annihilation 
operators as a p = u p a p + UpaiLp, where p ^ is the 
eigenvalue of the momentum operator p = 2nhxi/L, with 



n = (n x ,n y ,- 
the relations 



The real quantities u p and v p satisfy 



4 = 



V 



4) 



1/2 



2e E 



Ur,Vr> 



2e„ 



(1) 



and, consequently, the normalization condition u 2 — v 2 
1 . The energy of the elementary excitations entering ([ 
are given by the most famous Bogoliubov spectrum 
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where g = 47r7i 2 a/m is the coupling constant fixed by 
the s-wave scattering length a, and no is the condensate 
density which, at low T, coincides with the total density 
ri = N/V , apart from higher order effects in g arising 
from the quantum depletion of the condensate. At small 
p eq. (||) gives the phonon dispersion law e p = cp with 
the sound velocity given by c = yjgno/m. Notice that, 
due to the discretization of the values of p, the phonon 
dispersion law is ensured only if the condition Na/L ^> 1 
is satisfied. 

The total number of particles out of the condensate is 
given by the expectation value of the operator 



N n 



a l a p = H [ u p a p a p + v l a P a l 
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and takes the form 
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where 



f p = (a p a p ) = (expe p /k B T - 1) 
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is the number of elementary excitations present in the 
system at thermal equilibrium. Notice that the number 
of thermally excited atoms differs from the number of 
elementary excitations because of the factor (v,p + v 2 ). 
In the thermodynamic limit (V — ► oo) the sum (|4|) can 
be calculated by integration in momentum space. At low 
temperature this yields the result 



(N out ) = N 



m{k B T) 2 8^ 
\2ncfi 3 
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where N is the total number of particles. The term pro- 
portional to T 2 accounts for the phonon contribution to 
thermal fluctuations Jl(]] , while the temperature indepen- 
dent term gives the effect of quantum fluctuations 
The depletion of the condensate is obtained by simply 
taking the difference N = N— (N out ). For the ideal Bose 
gas confined in a box [u p = 1, v p = in (^J)] the same pro- 
cedure gives the well known result (N out ) — N(T/T C ) 3 / 2 , 
where k B T c = 2-nh 2 /m(n/C(3/2)) 2/3 is the Bose-Einstein 
transition temperature. 

In a similar way one can proceed to the calculation of 
the fluctuations of N out . Starting from (0) one finds 



(wL) = E [(« + v lf + Au >l) (fp + fp) 
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To calculate the expectation value of the terms with four 
quasi-particle operators appearing in (N 2 ut ) we have used 
Wick's theorem 



( a P a P a P ' a P' 



which holds for the creation and annihilation operators 
of statistically independent excitations. In the ideal Bose 
gas (u p = 1, v p = 0) result ((?]) is correct at all tempera- 
tures if one works in the grand canonical ensemble, and, 
as demonstrated in Ref. ||, holds also in the canonical 
ensemble except near and, of course, above T c . Physically 
this behaviour follows from the fact that the condensate 
plays the role of a "reservoir" for the particles out of the 
condensate. On the other hand, in interacting Bose sys- 
tems eq. (|7|) is valid only at low temperature where the 
number of excitations is small and their interaction can 
be neglected. Both in the interacting and non-interacting 
case eq. (0) allows us to calculate the fluctuations of the 
condensate in the canonical ensemble where one has triv- 
ially (<5iV 2 ) = (SN 2 out ) @. 

In equation (Q) one immediately points out the oc- 
currence of a strong infrared divergency arising from 
the term proportional to f 2 . In fact at low p one has 
u 2 ~ v 2 ~ gno/2cp and f p ~ ksT/cp. This results in a 
1/p 4 infrared divergency yielding a divergent integral in 
3D. Also in the ideal gas one finds a similar behaviour. In 
fact in this case u 2 = 1, v 2 = 0, while f p ~ 2mfcsT/p 2 . 



The occurrence of this divergency implies that the sum 
(0) cannot be calculated by integration in momentum 
space, but is dominated by the discretized sum over the 
low-energy phonon modes. The final result for the lead- 
ing terms in the large- N limit is given by 



(5N 2 ) = A 



mksT 



y4/3 



+ 



1 fc B Tm 2 c fVm 3 c 3 
73 Vl °S 



3tt 2 
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where A = 2/(vr) 4 £ n#0 1/n 4 = 0.105. The term con- 
taining log V arises from the contribution linear in f p of 
eq. (j?j) which yields an integral in momentum space ex- 
hibiting a logarithmic divergency. It is worth noting that 
the coefficient of the leading term does not depend on the 
velocity of sound (and hence on the interaction coupling 
constant), nor on the number of particles. In the ideal 
Bose gas (IBG) using cyclic boundary conditions the re- 
sult looks similar, but the coefficient differs by a factor 2 
11 
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In this case it is important to specify that the result holds 
in the canonical (C) ensemble. The IBG result (||) holds 
if the condition Na/L <C 1 is satisfied. 

The anomalous behaviour (||) exhibited by the fluc- 
tuations of the condensate reflects the occurrence of in- 
frared divergencies in the one-particle Green's function 
]l3"t . These phenomena are well accounted for by the hy- 
drodynamic theory of Bose superfluids (see for example 
Ref. |TJ]). Actually this theory can be directly employed 
to derive result (||) which consequently holds not only 
in the case of a weakly-interacting Bose gas, but also in 
strongly-interacting Bose fluids, like liquid 4 He. Notice 
that in this case the condensate fraction at T = siz- 
ably differs from the ideal gas value and that the sound 
velocity c, entering the log V term in eq. (||) , cannot be 
calculated using Bogoliubov theory. 

The fluctuations of the condensate can also be eval- 
uated at T — [fp — in eq. (0)]. In this case the 
sum has no infrared divergencies and can be safely cal- 
culated through direct integration in momentum space. 
This yields the result 



(5N 2 )(T = 0) = V^(an Q ) 3/2 V 



(10) 



revealing that, at T = 0, the fluctuations of the conden- 
sate have a "normal" behaviour, i.e. they are linear in 
V . These fluctuations are however large compared to the 
ones of the total number of particles, which in the grand 
canonical ensemble behave as (5N 2 ) ~ V 2 ^ 3 \ogV . This 
latter result is obtained by writing 



(SN 2 ) 



dr 1 dr 2 i'(\r2 - rJ 
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where v(r) is the density-density static correlation func- 
tion nv(\v2 — ri|) = (<5n(r 2 )5n(ri)). In an interacting 
Bose superfluid at T — the large r behaviour of v(r) is 
fixed by the asymptotic law |15fl 



for the fluctuations of the condensate in the canonical en- 
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Let us calculate the fluctuations of N in a sphere of radius 
R S> h/mc. The domain of integration in ( |lT| ) is then 
T\ ,f2 < R and one can write 



(SN 2 



Aim 
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dr r T{r)v{r) 
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where r(r) = 27ri? 3 (2 - 3r/2i? + r 3 /8i? 3 )/3 is the over- 
lapping volume of two spheres of radius R whose centers 
are separeted by the distance r. The integral of the r 
independent term in r(r) gives a contribution propor- 
tional to R 2 as can be seen by exploiting the sum rule 
4ir J Q dr r 2 v(r) = S(k = 0) = 0, holding for systems 
where the static form factor S(k) vanishes in the long- 
wavelength limit. The dominant contribution to (SN 2 ) 
comes thus from the term linear in r in r(r). One finally 
finds the relevant result fl^j 



(5N 2 )(T = 0) 



2nh 



-R*lo S (^) , (14) 
mc \ n ) 



holding with logarithmic accuracy. 

The analysis developed above can be easily extended to 
the case of a Bose gas confined in a harmonic potential of 



the form V Fr t = 



'/2, where, for simplicity, we have 



considered isotropic trapping. Also in this case one can 
define a proper "thermodynamic limit" . This is achieved 
by letting luq — > 0, N — > oo while keeping u>qN fixed 
p8| , ^9| . Differently from the gas confined in a box, in the 
presence of harmonic trapping the density of the system 
is inhomogeneous and it is convenient to use the following 
decomposition for the particle field operator 



^(r) = *(r) + + u< (r)aj) 



(15) 



where <£>(r) = (ip(r)) is the order parameter and the index 
i labels the elementary excitations of the system for which 
the normalization condition J dr(\ui\ 2 — |v;| 2 ) = 1 holds. 
The number of particles out of the condensate is given by 
the operator AT out = J dr(^(r) - $*(r))(i/'(r) - $(r)). 

Let us first discuss the non-interacting gas. In this 
case the fluctuations of N out can be written as (SN 2 ut ) = 
12ifi{fi + 1); where i labels the single-particle excited 
states of the oscillator Hamiltonian. By using the semi- 
classical approximation for these states p8| , we obtain 
the result [H 



6C(3) I T, 



(16) 



semble. Result (16) holds for temperatures smaller than 
the critical temperature for BEC [0 



k B T c = huj 
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which is kept fixed in the thermodynamic limit. The ap- 
plicability of the semiclassical approximation is ensured 
for temperatures much larger than the oscillator temper- 
ature (ksT 3> tiujo) and is a consequence of the fact 
that the sum J2i /»(/•» + 1) takes most of its contribution 
from excitations with energy ~ fcsT. The changes in 
the fluctuations of the condensate (|l6|) due to the use of 
the microcanonical ensemble can be also calculated and 
result in a change of the numerical prefactor . 

Result ( |l6| ) shows that the fluctuations of the conden- 
sate in a harmonic trap have a normal behaviour, i.e. 
they are proportional to N. This behaviour differs from 
result ([)]) for the ideal Bose gas in a box where the fluc- 
tuations have a stronger dependence on N (or V). This 
fact is not a surprise because it is well known that con- 
finement reduces the effects of thermal fluctuations. For 
example, the 2D ideal Bose gas in a harmonic trap ex- 
hibits Bose-Einstein condensation at finite temperature, 
differently from what happens in the uniform gas where 
thermal fluctuations have a destabilizing effect on the 
condensate. 

Interactions introduce however a deep difference in the 
above behaviour. In fact for repulsive forces the sys- 
tem exhibits a phonon-type behaviour at low excitation 
energies. This suggests that the fluctuations of the con- 
densate will be dominated by the discretized low energy 
modes which in the large- N limit obey the dispersion law 



e(n, I) = Hujo (2n 2 + 2n£ + 3n + () 
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These phonon-like collective excitations have been re- 
cently measured |2l]] and the observed frequencies are in 
excellent agreement with theory. Let us recall that the 
dispersion law ( |l8|) is valid if the conditions Na/auo ^ 1 
and e(n,£) <C n are satisfied. Here clho = y/fo/muJo is 
the oscillator length and [i is the chemical potential calcu- 
lated at zero temperature. The condition N a/ clho 3> lj 
which is well satisfied in the recent experiments on BEC, 
is the analogue of the condition Na/L ^> 1 which must 
be satisfied in the transition from the ideal to the Bo- 
goliubov gas in the box. Furthermore, the above con- 
dition yields the "Thomas-Fermi" result for the ground 
state density of the system no(r) = ^(r)! 2 = n/g(l — 
r 2 / R 2 )0{R c — r), where R c = 2^/(mw 2 ) is the radius of 
the condensate and /i = Hujo(l5Na/aHo) 2 ^ 5 /2. In the 
same limit one can also calculate the quantities itj and 
Vi of ( |l5| ) relative to the low- lying excitations (]l8|). One 
obtains the following leading behaviour |E2| 
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Ui(r) ~ -Vi(r) 




(19) 



where Xi{ r ) is the velocity potential associated with the 
collective mode satisfying the condition J drx* (r)Xj ( r ) = 
5ij. After a straightforward calculation one finds that 
the leading term in the fluctuations of the condensate 
exhibits an anomalous behaviour proportional to iV 4 / 3 : 



(20) 



where T c is given in (^) and B is a dimensionless coef- 
ficient 



{21+1) 



■f (2n 2 + 2n£ + Zn + £){2n n + 2n'£ + 3ri + I) 
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with 7 = 15 4/5 /(2(C(3)) 8/15 ). In the above equation 
n,n',£ = 0,1,2,.. and £ ^ if n(n') = 0, while 
ppn) _ p^( 2 ™) ^ are ^he polinomials introduced in Ref. 



It is worth noting that the inclusion of interaction 
reproduces a behaviour similar to the one exhibited by 
the gas confined in a box, i.e. proportional to iV 4 / 3 . This 
result reveals that the apparently similar behaviour in the 
ideal and interacting Bose gases in a box is accidental and 
that the physical effects induced by interactions show up 
in a clearer way in the presence of the harmonic trap. 

The consequences of interactions on the fluctuations of 
the condensate are expected to be even more dramatic in 
the case of attractive forces, where for values of N close to 
the critical size for the onset of instability, the monopolc 
vibration has a vanishing frequency with consequent ma- 
jor contribution to the sum (20). This situation, which is 
characterized by a full instability of the grand canonical 
ensemble, will be considered in a subsequent work. 

It is finally important to keep in mind that only the 
fluctuations of the condensate exhibit the anomalous be- 
haviour discussed in the present work. In particular the 
fluctuations of the density operator ip^(r)ip(r) are ex- 
pected to exhibit a normal behaviour in both uniform 
and non uniform interacting Bose gases. 
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